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ABSTRACT 

Aims. Understanding how galaxies evolve and assemble their mass across cosmic time is still a fundamental unsolved issue. To get 
insight into the various processes of galaxy mass assembly, the Mass Assembly Survey with SINFONI in VVDS (MASSIV) aims 
at probing the kinematical and chemical properties of a significant and representative sample of high-redshift (0.9 < z < 1.8) star- 
forming galaxies. 

Methods. This paper presents the selection function, the observing strategy and the global properties of the MASSIV sample. This 
sample contains 84 star-forming galaxies, selected from the VIMOS VLT Deep Survey (VVDS) and observed with the SINFONI 
integral-field spectrograph at the VLT. We present the redshift distribution, and derive the stellar masses and SED-based star formation 
rates (SFR). Integrated metallicities and the presence of type-2 AGNs are investigated using composite ID spectra built from VIMOS 
and SINFONI observations. 

Results. The MASSIV selection function, based on star formation criteria ([Oii]/13727 emission-line strength up to z ~ 1.5 and 
colors/UV absorption lines at higher redshifts), provides a good representation of "normal" star-forming galaxies with SED-based 
SFRs between 5 and 400 Mq yr"' in the stellar mass regime 10' - 10" M©. Analysis of typical emission-line ratios performed on 
composite spectra reveals that the contamination by type-2 AGNs is very low and that the integrated metallicity of the galaxies follows 
the well-known mass-metallicity relation. 

Conclusions. The MASSIV sample has been built upon a simple selection function, fully representative of the star-forming galaxy 
population at 0.9 < z < 1.8 for SFR > 5 Mq yr '. Together with the size of the sample, the spatially-resolved SINFONI data therefore 
enables us to discuss global, volume averaged, galaxy kinematic and chemical properties all accross the mass and SFR range of the 
survey to derive robust conclusions on galaxy mass assembly over cosmological timescales. All the data published in this paper are 
puWicly available at the time of publication following this link: http : //cosmosdb . lambrate . inaf . it/VVDS-SIMF OMI 
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1. Introduction 

The details of evolution processes driving galaxy assembly 
are still largely unconstrained. During the last decade, major 
spectroscopic and multi-wavelength photometric surveys (eg. 
VVDS, COSMOS, DEEP2, GOODS, CFHTLS, etc) explored 
in depth the high-redshift universe, allowing to follow the evolu- 
tion of large numbers of galaxies over cosmological timescales. 
Even if these data enabled significant progress concerning the 
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knowledge of the global properties (spectral energy distribution, 
colours, stellar mass and age, star formation rate, nebular metal- 
licity, size, etc) of different galaxy populations in different envi- 
ronments, more detailed constraints are needed to understand the 
main physical processes involved in the formation and evolution 
of galaxies. 

Based on observational and/or theoritical arguments, we 
know that galaxy merging (eg. Bournaud et al. 2007; Lin et al. 
2008; Conselice et al. 2008, 2009; de Ravel et al. 2009; Lopez- 
Sanjuan et al. 201 1) is at play in the growth of galaxies. Recently, 
cold gas accretion along cosmic filaments (Keres et al. 2005; 
Ocvirk et al. 2008; Genel et al. 2008, 2010; Dekel et al. 2009) 
has been emphasized as an efficient process to sustain high star 
formation rates and explain the clumpy nature of high-z disks. 
However observational signatures for cold accretion is still de- 
bated (e.g. Steidel et al. 2010; Le Tiran et al. 2011) and we do 
not know yet at which cosmic epoch one of these two processes 
(mergers and cold accretion) is dominant. 
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Spatially-resolved observations of high-redshift galaxies 
started to give some clues to answer this question as they allow 
to probe both the internal properties (kinematics, distribution 
of metals, ISM physics, etc) and close environment of distant 
galaxies. Powerful integral field spectrographs mounted of the 
largest ground-based telescopes (eg. GIRAFFE and SINFONI 
on the VLT, OSIRIS on Keck) recently made such studies possi- 
ble, giving access to the bright rest-frame optical emission lines 
such as Ua, H/3, [Nii]/t6584, [Oni]/I5007, and [Oii]/13727 in 
z ~ 0.4-4 star-forming galaxies. However, most of the previous 
surveys concentrated their efforts on the z > 2 universe (SINS: 
Forster Schreiber et al. 2006, 2009 - OSIRIS: Law et al. 2007, 
2009; Wright et al. 2007, 2009 - LSD/AMAZE: Maiolino et al. 
2008; Mannucci et al. 2009) with the exception of IMAGES who 
probed the kinematics of intermediate-z galaxies (eg. Puech et al. 
2008, 2010). However, much less is know in the z ~ 1-2 redshift 
range, except some very recent studies on small (< 30 galaxies) 
samples (Wisnioski et al. 2011; Mancini et al. 2011), where we 
know that L* galaxies, which are the major contributors to the 
cosmic star formation rate density at this epoch (Cucciatti et al. 
2011), are experiencing major transformations. Do cold gas ac- 
cretion be still an efficient process to assemble galaxies or do 
major mergers be predominant as it seems to be the case at later 
epochs (Yang et al. 2008, Keres et al. 2009)? 

MASSIV (Mass Assembly Survey with SINFONI in VVDS) 
tackles this issue by surveying a representative sample of star- 
forming galaxies in the redshift range z ~ 1-2. With the de- 
tailed information provided by SINFONI on individual galaxies, 
the key science goals of the MASSIV survey are to investigate 
in detail: (1) the nature of the dynamical support (rotation vs. 
dispersion) of high-z galaxies, (2) the respective role of mergers 
(minor and/or major) and gas accretion in galaxy mass assembly, 
(3) and the process of gas exchange (inflows/outflows) with the 
intergalactic medium through the derivation of metallicity gradi- 
ents. 

The MASSIV sample includes 84 star-forming galaxies 
drawn from the VIMOS VLT Deep Survey (VVDS) in the red- 
shift range 0.9 < z < 2.2. We stress the importance of working 
on representative and statistically-significant samples. Indeed, 
the acquisition of high-z galaxy kinematics is a long observa- 
tional process which requires to target galaxies that are repre- 
sentative of the main physical processes occurring at 1 < z < 2, 
in order to draw conclusions which are unbiased toward any 
specific population. We present the full MASSIV sample, fo- 
cusing on the selection criteria, the observing strategy, and the 
global properties derived mainly from SED fitting and integrated 
VIMOS and SINFONI composite spectra. The analysis of a first 
sample of 50 MASSIV galaxies is presented in related papers fo- 
cusing on the kinematical classification (Epinat et al. 2011), the 
evolution of scaling relations such as the baryonic TuUy-Fischer 
relation (Vergani et al. 2011), and the spatially-resolved metal- 
licity (Queyrel et al. 2011). The analysis of a sub-sample of 9 
MASSIV galaxies observed during pilot runs in 2005-2006 has 
already been published in Epinat et al. (2009) and Queyrel et al. 
(2009). The analysis of the full MASSIV sample, including the 
other 34 galaxies which have been observed with SINFONI, will 
be published in subsequent papers after the completion of data 
reduction and kinematical modeling. 

The paper is organized as follows. The selection of MASSIV 
targets is described in Section |2] The observing strategy with 
SINFONI is presented in Section |3] and includes a discussion 
on the advantages and drawbacks of using adaptive optics with 
the higher spatial sampling compared to seeing-limited obser- 
vations. In Section]?] we present the global properties (redshift. 



stellar mass, and star formation rate) of the MASSIV sample 
mainly derived from the SED fitting and we discuss how well it 
represents the z ~ 1-2 star-forming galaxy population. The inte- 
grated VIMOS and SINFONI composite spectra are shown and 
discussed in Section |5] The MASSIV sample global properties 
are compared to those of other Integral Field Unit (IFU) spec- 
troscopic samples at similar redshifts in Section |6] The paper is 
summarized in Section|7] 

Throughout this paper, we assume a standard A-CDM cos- 
mology, i.e. h - 0.7, Qm = 0.3 and Qa = 0.7 (Spergel et al. 
2003). For this cosmology, 1 "corresponds to ~ 8 kpc at z ~ 1 -2. 
Magnitudes are given in the AB-based photometric system, un- 
less explicitly stated otherwise. 



2. The MASSIV sample 

2.1. Parent VVDS catalogues and target selection criteria 

We have used the VIMOS VLT Deep Survey (VVDS) sample to 
select galaxies with known spectroscopic redshifts accross the 
peak of cosmic star formation activity (z ~ 1 - 2). 

The VVDS is composed of three /-band selected surveys to- 
talizing about 47 000 spectra of galaxies, quasars, and stars from 
z ~ to z - 5: (1) a Wide survey (17.5 < /ab < 22.5; Garilli 
et al. 2008), (2) a Deep survey (17.5 < /ab < 24.0; Le Fevre et 
al. 2004, 2005), and (3) an Ulti-a-Deep survey (23.00 < i'^^^ < 
24.75; Le Fevre et al., in prep.). Spectra have a spectroscopic 
identification at a confidence level higher than ~ 50%, 60%, 
81%, 97% and 99%, corresponding to the VVDS quality flags 
1,9, 2, 3 and 4. 

The VVDS is a complete magnitude-selected sample avoid- 
ing the biases linked to a priori color selection techniques. This 
sample offers the advantage of combining a robust selection 
function and secure spectroscopic redshifts. The latter are neces- 
sary to engage into long integration times on single galaxies with 
SINFONI being sure to observe the Ha line away from bright 
OH night-sky emission lines. From the existing VVDS dataset, 
we have access today to a unique sample of more than 4400 
galaxies in the redshift domain 0.9 < z < 2 with accurate and se- 
cure spectroscopic redshifts (Le Fevre et al. 2005). This sample 
being purely /-band limited, it contains both star-forming and 
passive galaxies distributed over a wide range of stellar masses, 
enabling us to easily define volume -limited sub-samples. 

For the MASSIV survey, we have defined a sample of 84 
VVDS star-forming galaxies at 0.9 < z < 2.2 suitable for 
SINFONI observations. Three selection criteria have been ap- 
plied successively. First, the MASSIV targets were selected 
to be star-forming galaxies, based on their [Oii]/i3727 or UV 
flux. The selection was performed on the measured inten- 
sity of [Oii]/l3727 emission line in the VIMOS spectrum (see 
Lamareille et al. 2009; Vergani et al. 2008) or, for the cases 
where the [Oii]/l3727 emission line was out of the VIMOS 
spectral range (i.e. for z > 1.5), on the UV flux based on 
their observed photometric UBVRIK spectral energy distribu- 
tion and/or UV rest-frame spectrum. The star formation crite- 
ria ensure that the brightest rest-frame optical emission lines 
(mainly Ha and [N ii]/l6584, or in a few cases [O iii]/15007) used 
to probe kinematics and chemical abundances, will be observed 
with SINFONI in the NIR J and H bands. Among these star- 
forming galaxy candidates, we have further restricted the sam- 
ple taking into account two important observational constraints: 
the observed wavelength of Ha line has to fall at least 9A away 
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from strong OH night-sky line^ in order to avoid heavy con- 
tamination of the galaxy spectrum by sky subtraction residuals; 
and the Ha line had to he in the J ox H band of SINFONI, 
thus prohibiting the redshift range z - 1.049 - 1.192. Finally, 
90% of MASSIV galaxies have been selected to be observable 
at higher spatial resolution with the adaptive optics (AO) system 
of SINFONI assisted with the Laser Guide Star (LGS) facility. 
In these cases, a bright star (/? < 18 mag) close enough to the 
target {d < 60") is needed for the zero-order tip-tilt corrections, 
as indicated in the SINFONI User Manual. 

In the following sub-sections, we detail the criteria used for 
the selection of MASSIV targets depending on the depth (Wide, 
Deep and Ultra-Deep) of the parent VVDS catalogue. 

2.1 .1 . The Wide VVDS sample 

The VVDS Wide used VIMOS at the ESO VLT to target four 
separate fields widely distributed on the sky (RA= 02h, lOh, 
14h, and 22h) and covering a total of 16 deg^. The spectro- 
scopic sample has been derived from an /-band selected photo- 
metric catalogue applying a pure apparent magnitude limit at /ab 
= 22.5. The final spectroscopic sample contains ~20000 galax- 
ies (broad-line QSOs are excluded on the basis of their peculiar 
spectral features) with a median redshift z ~ 0.58 (Garilli et al. 
2008). The average target sampling rates of VVDS Wide are be- 
tween 20% to 25% over 17.5 < /ab22.5, and are dependent of 
the magnitude (Garilli et al. 2008). The redshift success rates are 
~ 95, 70, 85, 80% for the 02h, lOh, 14h, 22h, respectively. 

MASSIV targets have been drawn from the VVDS-22h and 
VVDS-14h Wide fields. These fields contain a total of ~ 1500 
galaxies with a secure spectroscopic redshift over the redshift 
range 0.9 < z < 1.5. 

In order to select secure star-forming (SF) galaxies, we based 
our selection on the strength of [Oii]/l3727 emission lines mea- 
sured in VIMOS spectra up to z ~ 1 .5 (see Fig.[T]), when this line 
leaves the VIMOS spectral domain. The choice of the lower limit 
on [Oii]/l3727 equivalent width (EW) depends on the signal-to- 
noise ratio: it is equal to -25 A for S/N > 10 (this is the case for 
42 galaxies or ~ 66% of the full [Oii]^3727-selected MASSIV 
sample) and -40A for 6 < S/N < 10 (this is the case for 22 
galaxies). Galaxies with a S/N < 6 on [On]/13727 EW were not 
considered. This limit on [Oii]/l3727 EW ensures that the Ha 
emission line will be detected with SINFONI with a sufficient 
S/N ratio in a reasonable exposure time. A visual check has also 
been performed on VIMOS spectra to verify that the [Oii]/l3727 
line was not affected by a sky subtraction residual. We further ap- 
plied the OH sky lines and SINFONI bands criteria (see above). 

This selection strategy, based on [Oii]/l3727 equivalent 
width instead of using simply a flux limit, has been built on past 
experiences, especially on the detection rate obtained after the 
two pilot runs with SINFONI used to define the best selection 
strategy for the MASSIV Large Program. For the pilot runs, we 
used a selection based both on [Oii]/l3727 flux and equivalent 
width (see Fig. 1 in Epinat et al. 2009). After a careful inspec- 
tion of both detected and undetected targets with SINFONI, we 
concluded that a selection based on [Oii]/l3727 EW only allows 
to probe more "normal" star-forming galaxies by reaching SFR 
as low as a few Mq yr"' without decreasing significantly the de- 
tection rate. 



* this limit allows the detection of Ho- far enough from the galaxy 
dynamical center assuming rotating disks with a maximum circular ve- 
locity lower than 200 km s"' 



A total of 224 galaxies finally satisfied the above mentioned 
criteria in the VVDS-22h and -14h fields. In this sample, 183 
galaxies have a bright star close enough and could in principle 
be observed with AO-LGS. 

Finally, 21 galaxies have been selected in the VVDS Wide 
sample for SINFONI follow-up observations. The 13 MASSIV 
targets in the VVDS-22h field have been selected randomly 
among these 183 galaxies, but 4 only were observed with 
the AO-LGS facility. The 8 additionnal targets selected in the 
VVDS-14h field do not have any bright star close enough to be 
observed with AO-LGS. 

2.1 .2. The Deep VVDS sample 

The VVDS Deep survey used VIMOS at the ESO VLT to 
target data in the VVDS-0226-04 (Le Fevre et al. 2005) and 
VVDS-CDFS fields (Le Fevre et al. 2004). In this study, we 
use the VVDS-0226-04 field only, which covers 2200 arcmin^ 
of sky area. It includes 9842 spectra described in Le Fevre et 
al. (2005), plus 2826 spectra acquired later with the same set-up 
. Observations of several z > 1.4 targets, with assigned VVDS 
quality flags of 0, I, or 2, have been reapeated in order to as- 
sess their real redshift distributions as detailled in Le Fevre 
et al. (in prep). The measurement of the target sampling and 
redshift success rates are estimated in Ilbert et al. (2005) and 
Cucciati et al. (2011), the later work is accounting for the re- 
peated observations. In summary, the average target sampling 
rate is ~ 24 - 25%, and is calculated as a function of the projec- 
tion of the angular size of each target in the .jc-axis of the image. 
The redshift success rate is estimated as a function of /ab and 
redshift; for instance at z < 1.5, it is ~ 100% at /ab < 22.5%, 
and it smoothly decreases down to ~ 80% at the faintest mag- 
nitudes. Stellar masses have been derived for all galaxies with a 
measured redshift (see eg. Walcher et al. 2008). The VVDS flux 
limit translates at different redshifts into dififerent lower luminos- 
ity limits. This, in turn, translates into a broad mass selection cut 
at each redshift, reflecting the scatter in the Mass-Luminosity 
relation. Meneux et al. (2008) have studied the stellar mass in- 
completeness as a function of redshift. At z < 1.5, it is on aver- 
age 84% complete in stellar mass for M»> IO^-^Mq. 

MASSIV targets have been drawn from this Deep VVDS 
field which contains a total of ~ 2600 galaxies with secure spec- 
troscopic redshift over the redshift range 0.9 < z < 2.0. 

The criteria applied for the selection of secure star-forming 
galaxies is identical to the one used for the Wide VVDS cata- 
logue, i.e. based on the strength of the [O n]/l3727 emission line. 
We further applied the OH sky lines and SINFONI bands criteria 
as for the Wide fields. A total of 490 galaxies satisfy these cri- 
teria in the VVDS-02h Deep field. In this sample, 300 galaxies 
have a bright star close enough and could in principle be ob- 
served with AO-LGS. 

Finally, 29 MASSIV targets have been selected randomly in 
this parent sample of 300 galaxies. Among these targets one is 
observed with the AO-LGS facility. 

2.1 .3. The Ultra-Deep VVDS sample 

The VVDS Ultra-Deep has assembled a sample of 1200 new tar- 
gets over 576 arcmin^ within the 2200 arcmin^ sampled by the 
Deep survey. Thanks to its large spectral range (~ 3600-9350A) 
and long exposure time (~ 15 hours), this unique spectroscopic 
survey enables one to securely identify high-redshift galaxies up 
to z ~ 5, filling the so-called "redshift desert" between z = 1 and 
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Fig.l. Selection of MASSIV 0.9 < z < 1.5 stai--forming 
galaxies with secure redshift and measured [On]/i3727 emis- 
sion line in the VVDS Wide (VVDS-14h and -22h, Iab < 
22.5, small triangles) and Deep/Ultra-Deep VVDS-02h (/ab 
< 24/24.75, small squai-es) fields for SINFONI follow-up ob- 
servations. Among these galaxies we selected those showing 
strong enough [Oii]/l3727 emission line for Ha to be easily de- 
tected in the NIR, and satisfying the OH and AO/LGS criteria 
(filled symbols, see text for details). The 63 galaxies selected 
for the MASSIV survey in the redshift range 0.9 < z < 1.5 are 
indicated as blue triangles (VVDS-14h and -22h Wide fields), 
red squares (VVDS-02h Deep field), and green circles (VVDS- 
02h Ultra-Deep field). The distribution of [Oii]/l3727 equivalent 
width (EW) and flux are projected onto the x- and y-axis re- 
spectively. These histograms are normalized to the total num- 
ber of objects. Empty-dot histograms are for the full VVDS 
Wide, Deep, and Ultra-Deep samples with measured [Oii]/l3727 
emission line (A^ = 2397, median EW = -38A, median flux 
= 4.4 X 10^'^ergs s"' cm"^). Empty-solid histograms are for 
the VVDS Wide, Deep, and Ultra-Deep sample galaxies com- 
plying with the selection criteria ([Oii]/l3727 strength, OH sky 
lines and redshift range, see text for details, = 728, median 
EW = -54A, median flux = 6.9 x 10^'^ergs s^' cm^^) pjjjg^j 
cyan histograms are for the selected MASSIV galaxies (A^ = 63, 
median EW = -61A, median flux = 8.3 x 10^'^ergs s"' cm"^). 
The dashed line indicates the selection box of MASSIV targets 
based on [Oii]/l3727 EW, flux and spectrum SNR (see text for 
details) The cross in the upper-left corner shows the typical error 
bai-s on [O n]^3727 flux and EW. 



2. The average target sampling rate of the Ultra-Deep survey is 
constant, i.e. ~ 4% at 23 < Iab < 24.75. Combining the VVDS 
Deep and Ultra-Deep surveys (17.5 < /^g < 24.75), Cucciati 
et al. (2011) have derived the following rates. The target sam- 
pling rate starts from ~ 20% at Iab ~ 22.5, reaches ~ 30% at 
Iab ^ 23, and then decreases to 2% at Iab ^ 24.75. The redshift 
success rate is measured as a function of magnitude and redshift. 
From z = 1.5 to z = 2, it decreases from ~ 80% to ~ 65% (see 
Le Fevre et al., in prep.). 



This Ultra-Deep VVDS catalogue has been used to select 
MASSIV galaxies in the redshift domain z ~ 1.24 - 1.8, in order 
to target Ho' in the H band. It complements nicely the ~ 1 .2 - 1 .5 
redshift range, already populated with sources drawn from the 
VVDS Wide and Deep, and allows us to extend the MASSIV 
sample up to z ~ 1.8 with a robust selection. 

The selection procedure to identify secure star-forming 
galaxies for the SINFONI follow-up is different from the one de- 
scribed above for the Wide and Deep VVDS catalogues. Indeed, 
when the MASSIV target selection has been performed, no sys- 
tematic measurement of spectral features (emission and absorp- 
tion lines) was applied to the Ultra-Deep sample. We thus se- 
lected galaxies with either a prominent [O ii]/l3727 emission line 
(as measured with IRAF splot task) or, for redshifts greater than 
~ 1 .5, a rest-frame UV spectrum typical of star-forming galaxies 
(i.e. with strong Civ/11550 absorption line and/or blue contin- 
uum, see Fig.|7]). We further applied the OH sky lines criterium 
as for the Wide and Deep fields restricting the VVDS Ultra-Deep 
sample to 86 galaxies fulfilling our criteria. In this sample, 55 
galaxies have a bright star close enough and could be in princi- 
ple observed with AO-LGS. 

Finafly, 34 MASSIV targets (including 18 [On]^3727- 
selected at z < 1.5) have been drawn from this parent sample 
of 55 galaxies. Among these targets 6 are observed with the AO- 
LGS facility. 

2.1 .4. The final MASSIV sample 

The full MASSIV sample is extracted from the three VVDS 
samples described above, for which the completeness func- 
tions are well-controlled. This sample, selected either on the 
[On]/i3727 EW strength or on the UV intensity, is essentially 
limited in flux translating into a lower limit for the star forma- 
tion rate (see Fig. [6]). Figure [T] illustrates the selection proce- 
dure based on [On]/l3727 emission line strength for the Wide, 
Deep, and Ultra-Deep VVDS samples. The distribution of the 63 
MASSIV targets, selected among the strongest [O ii]/l3727 emit- 
ters (lower limit on EW only) span a wide range of [Oii]/l3727 
flux (from ~ 1.2 to 45 X 10"'^ ergs s"' cm"^) and EW (from 
~ -200 to -30A). The median values of these distributions 
(EW = -61A, flux = 8.3 X 10^'^ergs s"' cm^^) sHghtly 
higher than the median values of the combined VVDS Wide, 
Deep and Ultra-Deep parent samples (EW - -54A, flux = 
6.9 X 10^'^ergs s"' cm"^) complying with the selection crite- 
ria listed above ([On]/l3727 strength and redshift range). It is 
important to stress that the expected star formation rates (rang- 
ing from ~ 8 to ~ 1000 M© yr"' with a median value of 41 
Mq yr"', and using the [Oil] -based calibration of Argence & 
Lamareille 2009) corresponding to these limits imposed on the 
observed [Oii]/l3727 line strength are in very good agreement, 
taking the dust attenuation into account, with the SED-based star 
formation rates derived in Section|42] 

The properties of the MASSIV sample are further described 
in the next sections focusing on the selection function (sect. |4.3| l 
and the comparison with other major IFU surveys of high-z 
galaxies (sect.|6]). 

3. SINFONI observations 

3.1. Strategy 

The observations of the MASSIV sample were performed with 
SINFONI (Eisenhauer et al. 2003; Bonnet et al. 2004) mounted 
at the Cassegrain focus of the VLT UT4 telescope. Most of the 
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data (90%) have been collected in service mode during several 
observing periods between April 2007 and January 20 11, as part 
of the ESO Large Program 179-A.0823. However, eleven galax- 
ies of the MASSIV sample (including three with duplicated ob- 
servations) were observed during two pilot runs with SINFONI 
(September 2005 and November 2006, in visitor mode). The 
analysis of nine of these galaxies has been published in Epinat 
et al. (2009) and Queyrel et al. (2009). 

To map the Ha, [Nii]^6584, and [S ii]/l/l6717,6731 emission 
lines, or [Oiii]/I5007 and H/3 for four galaxies of the MASSIV 
sample, we used the J or H gratings, depending on the redshift of 
the sources, which have a nominal FWHM spectral resolutions 
ofR~ 1900 and 2900 respectivel}]^ The vast majority of the ob- 
servations (74 galaxies) were carried out in seeing-limited mode 
with the largest pixel scale of 0.125"/pixel giving a Field-of- 
View (FoV) of 8"x8". We observed a total of eleven targets with 
adaptive optics (AO) assisted with the Laser Guide Star (LGS). 
For seven of them, we selected the intermediate 0.05"/pixel scale 
with FoV of 3.2"x3.2"to take full advantage of the gain in an- 
gular resolution provided by the AO, achieving a FWHM resolu- 
tion of ~ 0.25" . For the four other targets with AO data, we used 
the larger pixel scale as trade-off between enhanced angular res- 
olution and sensitivity (see below for a detailled discussion). The 
observing conditions were generally good, with clear to photo- 
metric sky transparency and typical seeing at NIR wavelengths 
with FWHM = 0.5 -0.8". 

Depending on the FoV, we adopted two different observing 
strategies to maximize the efficiency of observations. When the 
largest 8"x8" FoV was used, i.e. for the vast majority of ob- 
servations, we applied an "object nodding" technique by mov- 
ing the target from one comer to the other of the FoV. For 
the seven galaxies observed with AO-LGS at 0.05"/pixel, we 
adopted the strategy of taking the sky frames away from the tar- 
get because the source size did not allow nodding within the 
small 3.2"x3.2"FoV. In such cases, the telescope was pointing 
between the object (O) and adjacent sky regions (S) in an O-O- 
S-O-O-O-O-S-O-O pattern. 

The individual exposure times vary between 300s and 600s 
depending mainly on the grating J or H used for observations 
(individual exposures of 900s have been used for a few galaxies 
during the pilot runs, see Epinat et al. 2009). The total on-source 
integration times range on average from 80 to 120 minutes. The 
total integration times were driven by the [O ii]/l3727 emission 
line intensity as measured in VIMOS spectra or, when this infor- 
mation was not available, by the /-band magnitude of the target. 

In order to estimate the spatial PSF, we obtained exposures 
of a bright star close to the target (or tip-tilt stars for AO-LGS), 
used also to center the galaxy into the FoV with a blind offset. 
For flux calibration and atmospheric transmission correction, we 
observed telluric standard stars as part of the standard ESO cal- 
ibration plan. For details on the data reduction procedures, we 
refer to Epinat et al. (201 1). 

3.2. Consistency between AO and seeing-limited 
observations: VVDS220596913 as a test case 

As stated previously, even if the majority of MASSIV galaxies 
have been observed with SINFONI in natural seeing conditions, 
eleven galaxies of the sample benefit from AO-LGS observations 
in order to achieve a better spatial resolution. However, these AO 



^ The nominal spectral resolution is slightly lower with the 
0.05"/pixel scale used for some AO-assisted observations (see p. 68 of 
the SINFONI User Manual, version P89). 



observations were not performed always with the same spatial 
sampling. Indeed, the observing strategy we adopted with the 
small pixel scale (0.05") at the beginning of the survey requires 
the acquisition of sky frames in separate exposures (see above). 
Thus, for a given total observing time, only 80% can be spent 
on the object. To optimize the survey efficiency, we decided to 
modify our strategy for AO-assisted observations from Period 
82, shifting to a larger pixel scale (from 0.05" to 0.125"). Thus, 
at the end, seven galaxies out of eleven have been observed with 
the smallest pixel scale. For an homogeneous analysis and in- 
terpretation of the resolved properties of MASSIV galaxies, it is 
then important to check whether the differences in the final spa- 
tial resolution/sampling achieved with the SINFONI datacubes 
have an impact on the derived galaxy properties (Ha flux, veloc- 
ity fields, etc). 

The galaxy VVDS220596913 at z ~ 1.27 is a perfect 
"test case" for this purpose as it has been observed twice with 
SINFONI. A natural seeing observation has been first acquired 
in 2007 as part of a pilot run (see Epinat et al. 2009). This 
object was then re-observed in the frame of the MASSIV sur- 
vey, but this time using AO-LGS. Data from the pilot run have 
been re-reduced using the same method than for the MASSIV 
data (Epinat et al. 2011). In Table [T] we summarize the ob- 
serving conditions and derived parameters for the seeing-limited 
and AO-assisted SINFONI observations of VVDS220596913. 
Taking into account both the variation in pixel surface and the 
difference in exposure time between the two observations, we 
expect the signal-to-noise per pixel of the AO observation to be 
lower than in the seeing-limited one by a factor of 2.9. Using 
the Ha flux maps shown in Figure [2] we measure a factor of 
~ 2.5. This slight difference may be due to the fact that the sig- 
nal is less diluted in AO observations and/or observations were 
not performed under the same photometric conditions. 

From our AO observations with high spatial sampling 
(0.05"/pixel), it was possible to resample the data in order to in- 
fer what would be the gain in using AO with a larger pixel scale 
with respect to non-AO observations and with respect to small- 
pixel AO observations. In order to avoid any interpolation, we 
used a spatial resampling of 0.10" instead of 0.125". We also 
mimic seeing-limited observations by adding a blurring of the 
datacube to reach the same final spatial resolution of our non-AO 
observation but with a 0.10"/pixel spatial sampling. We created 
datacubes from which we derived Ha flux and kinematics maps 
after a two pixels FWHM spatial Gaussian smoothing, as is the 
case for the observed data. The limitation of these simulations is 
that the lowest surface brightness regions could be not detected 
at all, thus, rebinning the data would not be strictly equivalent to 
real observations. This is further supported by the fact that real 
and simulated non-AO observations are not strictly identical in 
the relative flux of the various components (see Figure|2]i. 

The gain of AO-assisted observations with the large pixel 
scale with respect to non-AO observations is clearly visible. 
Indeed, in the non-AO observation it is difficult to distinguish the 
four blobs/clumps that are clearly separated in AO-assisted ob- 
servations independently of their pixel scale (see Figure|2]i. Due 
to surface brightness detection limit, these clumps might be even 
more clearly detected with the large pixel scale (both fainter sur- 
face brightness limit and longer exposure time). However, the 
interpretation would still remain the same since we do not de- 
tect significant change in the velocity field which remains mono- 
tonic. 

We further compared the size of the Ha emission, the flux in 
this line as well as, assuming a rotating disk model (see Epinat 
et al. 2011), the kinematics position angle of the major axis. 
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Fig. 2. Ha maps and velocity fields of VVDS220596913 (z ~ 1.27) obtained with and without adaptive optics. From top to bottom: 
Ha flux map, velocity field and velocity dispersion map (corrected for spectral resolution). From left to right: AO-assisted SINFONI 
observation (pixel size = 0.05"), simulated observation with AO and a pixel size of 0. 10", simulated seeing-limited observation with 
a pixel size of 0.10", seeing-limited SINFONI observation (pixel size - 0.125"). The circle in the lower left corner of each image 
represents the effective FWHM spatial resolution (a two pixel Gaussian smoothing has been applied to each cube before map 
extraction). The color scale of the Ha flux maps is arbitrary. A threshold of S/N per pixel > 3 is used for the display of velocity 
maps. North is up and East is left. 

Table 1. VVDS220596913: physical parameters derived from SINFONI observations (with and without AO) and simulations. 
Parameter cr corresponds to the intrinsic velocity dispersion corrected for both instrumental broadening (spectral resolution) and 
beam smearing (see Epinat et al. 201 1 for details) 



Parameter 


SINFONI 


Model 


Model 


SINFONI 




AO 


Binned 


Binned and blurred 


Seeing-limited 


Spatial sampling ["] 


0.05 


0.10 


0.10 


0.125 


Spatial resolution ["] 


0.183 


0.183 


0.573 


0.573 


Exposure time [minutes] 


120 






160 


Ha flux [xlO""' ergs s"' cm"^] 


3.0 






4.0 


Ho- size ["] 


2.6 


2.5 


2.7 


3.4 


Position angle [°] 


247 ± 1 


250 ± 1 


248 ± 1 


246 ± 1 


V„,„.j [km s"'] 


141 + 2 


127 ±2 


155 ±2 


165 ±3 


cr [km s '] 


38 ±28 


41 ±24 


37 ±22 


66± 22 



the maximum rotation velocity and the mean velocity disper- 
sion (see Table [T]i. We find that all these quantities are lower in 
the AO observation. However, except for the velocity dispersion, 
they remain compatible within 25%. This can be understood by 
the fact that the AO observation is less deep and thus misses 
some low-surface brightness regions. Using a similar signal-to- 
noise threshold, we end up with final maps that are slightly less 



extended than in the seeing-limited observations. Note that we 
also find a very good agreement between the various kinemat- 
ics position angles of the major axis obtained from kinematics 
modeling. 

We have perfomed advanced tests to better understand the 
larger variations observed in the velocity dispersion. First, it has 
to be noticed that the difference before correction for spectral 
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PSF is only of the order of 10 km s . The accuracy of the 
wavelength calibration could thus be a cause for the broaden- 
ing of lines in non-AO observations due to the nodding strategy. 
However, we have checked that the wavelength calibration can- 
not account for the observed difference (it could only account 
for ~ 1 km s"'). The determination of the spectral PSF is also 
unlikely to explain alone the difference. One other explanation 
could be a spatial offset between various individual exposures 
since several components at various velocities can be summed 
up if spatial offsets are not correctly controlled. However, we 
would expect the same offsets in AO and non-AO observations. 
The accuracy of the offsets cannot be checked since no bright 
component is present in the data. Two other effects could explain 
the observed discrepancy in velocity dispersions: the lower SNR 
in the AO data and the presence of a sky line residual close to 
the Ha line. The SNR has been probed using a single OB for the 
non-AO data, but is seems that the velocity dispersion remains 
at the same value, even if locally it can change by about 5 to 10 
km s ' . However it could be that low surface brightness regions 
are not observed in the AO observation. The final difference is 
probably a combination of all these effects. It has however to be 
mentionned that the final velocity dispersion remains compatible 
between the different configurations within the computed uncer- 
tainties. 

Based on this test case, we conclude that the final resolu- 
tion/sampling of SINFONI observations does not impact signif- 
icantly the derived Ha flux and velocity maps. 

4. Global properties of the MASSIV sample 

4.1. Redshift distribution 

The redshift distributions of the MASSIV sample are compared 
with those of the parent VVDS catalogues in Figu re [3] Median 
values and associated dispersion are listed in Table PTMASSIV 
galaxies are distributed in the redshift range between z ~ 0.94 
and 1.80, with the exception of VVDS220148046 which has a 
redshift z = 2.24. This galaxy was supposed to be in the same 
redshift range as other MASSIV targets but it appeared that the 
redshift derived from the VIMOS spectrum was wrong. The red- 
shift of this galaxy is now secured with the measurement of H/S 
and [Oiii]/I5007 emission lines in the SINFONI spectra. The 
median redshift for the MASSIV sample is z = 1.33. The low- 
est redshift bin (z ~ 1) is populated with galaxies selected in 
the Wide and Deep VVDS samples, whereas most of the galax- 
ies at higher redshifts (z > 1.5) are selected from the Ultra- 
Deep VVDS sample only. The median redshift values for the 
MASSIV galaxies drawn from the Wide, Deep and Ultra-Deep 
VVDS samples are z = 1 .05, z = 1 .28, and z- 1 .48 respectively. 

4.2. Physical parameters derived from SED fitting 

In addition to optical VIMOS spectroscopy, all VVDS fields 
have a multi-wavelength photometric coverage: B, V,R,I from 
the VIRMOS Deep Imaging Survey (VDIS, Le Fevre et 
al. 2004), M* g', r', i', z from the Canada-France-Hawaii 
Telescope Legacy Survey (CFHTLS, but for VVDS-22h and - 
02h only), 7, H, and K from UKIDSS, and a partial coverage 
in J and K bands from VDIS (lovino et al. 2005). The VVDS 
Deep field has also been covered in the IRAC 3.6, 4.5, 5.8, and 
8.0 jum bands by Spitzer within the SWIRE survey (Lonsdale et 
al. 2003). 

Making use of the optical spectra and of all the photom- 
etry available, we have fitted the spectral energy distribution 




Redshift 

Fig. 3. Redshift distribution of the MASSIV samples (filled his- 
tograms) compared with the parent (i.e. restricted to the 0.9-2.0 
redshift range) VVDS samples (white histograms). From top 
to bottom: complete (magenta). Wide (green). Deep (red) and 
Ultra-Deep (blue) MASSIV and VVDS samples. Histograms are 
arbitrarily normalized. The redshift distributions of MASSIV 
reflect the selection criteria, but are also importantly affected 
by the observability of the target emission lines (mainly Ha, 
[Oiii]/I5007 in a few cases) in the NIR atmospheric bands 
(explaining the gap around z ~ 1.1) and between the OH 
night sky lines. MASSIV galaxies are distributed in the red- 
shift range between z ~ 0.94 and 1.80, with the exception of 
VVDS220148046 which has a redshift z = 2.24 (see text for 
details). 



(SED) of each galaxy to obtain estimates of the stellar mass 
and star formation rate (SFR) of our sample and of the parent 
VVDS sample. The fitting of the SED has been performed us- 
ing the GOSSIP spectral energy distribution modeling software 
(Franzetti et al. 2008): photometric and optical spectroscopic 
data were fitted with a grid of stellar population models, gen- 
erated using the BC03 population synthesis code (Bruzual & 
Chariot 2003), assuming a set of "delayed" star formation histo- 
ries (see Gavazzi et al. 2002 for details), a Salpeter (1955) IMF 
with lower and upper mass cutoffs of respectively 0.1 and 100 
Mq, a metallicity ranging from 0.02 and 2.5 solar metallicity 
and a Large Magellanic Cloud reddening law (Pei 1992) with 
an extinction E{B - V) ranging from to 0.3. The parameters 
for the best-fitting model for each galaxy are taken as the best 
fitting values for both the galaxy stellar mass and SFR. On top 
of the best fitting values, GOSSIP computes also the Probability 
Distribution Function (PDF), following the method described in 
Walcher et al. (2008). The median of the PDF and its confidence 
regions are then used to derive a robust estimate of the parameter 
and associated uncertainty that is to be determined. 

We have used all the photometric information available for 
each galaxy within the MASSIV sample. Each fit has been vi- 
sually inspected, and in few cases photometric points strongly 
deviating from the global fit have been removed. Within VVDS, 
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Table 2. Statistics on the redshift and global physical parameters of the MASSIV and parent VVDS samples 



Sample 




Redshift 




log(Stellar Mass) [Mq] 




log(SFRsED) [Moyr-'] 






MASSIV 




VVDS 




MASSIV 


VVDS 




MASSIV 


VVDS 






Median ± a 


A' 


Median ± a 


A' 


Median ± cr 


A' Median ± cr 


N 


Median ± cr N 


Median ± cr 


A' 


Total 


1.33 ±0.13 


84 


1.06 ±0.13 


4446 


10.15 ±0.30 


84 10.51 ±0.43 


4248 


1.50 ±0.31 84 


1.22 ±0.35 


3023 


Wide 


1.05 ±0.14 


21 


I.00±0.09 


1483 


10.68 ± 0.24 


21 10.89 ±0.34 


1409 


1.80 ±0.25 21 


1.36 ±0.27 


783 


Deep 


1.28 ±0.15 


29 


1.09 ±0.13 


2643 


10.08 ±0.14 


29 10.30 ±0.38 


2537 


1.37 ±0.28 29 


1.12 ±0.36 


1982 


Ultra-Deep 


1.48±0.I1 


34 


1.33 ±0.25 


320 


10.15 ±0.29 


34 10.21 ±0.33 


302 


1.50 ±0.27 34 


1.33 ±0.47 


258 




9 10 11 12 

log(Slellar Mass) [Msun] 



Fig. 4. Stellar mass distribution of the MASSIV samples (filled 
histograms) compared with the parent VVDS samples (white 
histograms). From top to bottom: complete (magenta). Wide 
(green), Deep (red) and Ultra-Deep (blue) MASSIV and VVDS 
samples. Histograms are normahzed to the total number of ob- 
jects. 




log(SFR,,„) [Msun/yr] 

Fig. 5. SED-derived SFR distribution of the MASSIV samples 
(filled histograms) compared with the parent VVDS samples 
(white histograms). From top to bottom: complete (magenta). 
Wide (green). Deep (red) and Ultra-Deep (blue) MASSIV and 
VVDS samples. Histograms are normahzed to the total number 
of objects. 



spectroscopic data are flux calibrated at the 10-20% level (Le 
Fevre et al. 2005). In order to fit them together with photome- 
try, we have renormalized spectroscopic data to the photometric 
;'-band magnitude. As MASSIV galaxies have been selected on 
the basis of their star-forming activity, we have been able to fur- 
ther constrain the fitting by imposing that the resulting model 
should not have a Z)4ooo break larger than 1 .25 with some tol- 
erances to account for the uncertainties on the measurement of 
^^4000 break. This is a reasonable value that excludes early-type 
passive galaxies but keeps the possibility to deal with an early- 
type star-forming galaxy (see e.g. Vergani et al. 2008). 

We have applied the same procedure also to the VVDS par- 
ent samples described in Section 2.1 all available photometry 



has been used for each galaxy, together with spectroscopy, renor- 
malized to the / '-band magnitude, or to the /-band magntitude 
when the previous one was not available (as it is the case for the 
VVDS-14h field). As the parent sample comprises all kinds of 
galaxies, we have not imposed a piior on the D4000 break. 

Stellar mass and star formation rate derived from the SED 
fitting (median values of the PDF and associated uncertainties) 
for MASSIV galaxies are listed in Table|5] 



The stellar mass distributions of the MASSIV sample are 
shown in Figure |4] Median values and associated dispersion are 
listed in Table |2] MASSIV galaxies are distributed in the stellar 
mass range between ~ 3 x 10^ and 6 x 10'' M©. The median 
stellar mass for the MASSIV sample is 1.4 x 10'" Mq. The me- 
dian stellar masses for the MASSIV galaxies drawn from the 
Wide, Deep and Ultra-Deep VVDS samples are 4.8 x 10'° Mq, 
1.2 X 10'° Mq, and 1.4 x 10'° Mq respectively. 

The SED-derived star formation rate distiibutions of the 
MASSIV sample are shown in Figure l5l Median values and as- 
sociated dispersion are listed in TableHl MASSIV galaxies are 
distributed in the SFR range between ~ 5 to 400 Mq yr"'. The 
median SED-derived SFR for the MASSIV sample is ~ 32 Mq 
yr"'. The median star formation rates for the MASSIV galaxies 
drawn from the Wide, Deep and Ultra-Deep VVDS samples are 
~ 63, 23, and 31 M© yr"' respectively. 

The relation between the SED-based SFR and stellar mass is 
shown in Figure |6] The location of MASSIV galaxies selected 
in the Wide, Deep and Ultra-Deep catalogues is compared to the 
one of the parent VVDS sample. Most of the MASSIV galaxies 
occupy an area bounded by the star formation "main sequences" 
defined empirically at z ~ 1 and z ~ 2 (Bouche et al. 2010). This 
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Fig. 6. SED-derived star formation rate as a function of stel- 
lar mass. The MASSIV sample (Wide - green circles, Deep - 
red squares and Ultra-Deep - blue triangles) is compared to the 
VVDS parent sample (small filled squares). The lines represent 
the empirical relations between SFR and stellar mass for differ- 
ent redshifts between z = and z = 3 following the analytical 
expression given in Bouche et al. (2010). 



figure shows cleary that the MASSIV sample is representative of 
the overall population of star-forming galaxies at z ~ 1-2. 

4.3. Discussion on seiection effects 

Having been assembled using different flux-limited VVDS cat- 
alogues, it is worth assessing what part of the z ~ I - 2 pop- 
ulation is represented by our MASSIV sample with respect to 
the unbiased population of high-redshift galaxies probed by the 
VVDS. Indeed, VVDS is a purely apparent magnitude-selected 
sample, without any color selection unlike color-selected sam- 
ples produced in the same redshift range as MASSIV (e.g. BzK, 
BM/BX). The VVDS sample is thus representative of the overall 
star-forming population of galaxies at high redshifts. 

The most stringent selection criteria used to built the 
MASSIV sample is certainly the requirement of a minimum 
[Oii]/l3727 equivalent width (see Sect, 
limit is likely to translate into an overall ' 



2.1 



This sensitivity 
"bias" towards younger 
and more actively star-forming systems. To check for this possi- 
ble effect, we compare in Figure[T|the [Oii]/l3727 flux and EW 
distributions of the final MASSIV sample to those derived for 
the VVDS parent catalogues. We clearly see that the selection 
criteria based on the [Oii]/i3727 EW ends up with a sample of 
parent VVDS galaxies (empty-solid histograms in Fig. [T} with 
higher median values of the [Oii]/l3727 EW (-54A compared 
with -38A) and flux (6.9 compai-ed with 4.4 x 10"'^ ergs s"' 
cm"^). However, the final random selection of MASSIV targets 
(filled cyan histograms in Fig.[T]) among this parent sample does 
not introduce any additional strong selection effect as the median 
values of [On]/13727 EW and flux (-61 A and 8.3 x lO"'^ ergs 
s~' cm"^) do not vary significantly, even if we slightly tend to se- 



lect the strongest [On]/l3727 emitters. A Kolmogorov-Smirnov 
(K-S) statistical test performed on the [On]/i3727 flux distribu- 
tions of the MASSIV and [Oil] -selected VVDS samples further 
shows a high probability for these two samples to be drawn from 
the same distribution. 

With a selection based on [On]/13727 emission-line or UV 
intensity, MASSIV is thus essentially a flux-limited sample, 
translating into a SFR lower limit for a given stellar mass (see 
Fig.[6|l, and reaching SFRs as low as a few Mg yr '. 

Looking at the stellar mass distributions in Figure]?] and me- 
dian values listed in Table ]2] we see clearly that we are prob- 
ing different mass regimes depending on the considered VVDS 
parent sample. The Wide VVDS sample, with its limiting mag- 
nitude /ab ^ 22.5, contains more massive galaxies than the 
Deep/Ultra-Deep samples (/ab < 24 and 24.75 respectively). 
This is due both to the larger area and shallower limiting mag- 
nitude of the wide survey compared with the deeper ones. But 
the important point is that the MASSIV sample, drawn from the 
Wide, Deep, and Ultra-Deep VVDS samples, is probing a large 
and representative range of stellar masses allowing us, in fine, to 
draw robust conclusions on the nature of the overall population 
at high redshifts. The part of the VVDS z ~ 1-2 population that 
is most clearly absent among the MASSIV sample is the mas- 
sive quiescent tail at low star formation rates (see Figures ]4] and 
]5]l. Indeed such objects would be difficult to detect as no or very 
faint [On]/i3727 is expected, at least from star formation. 

We conclude from this section that in spite of the selection 
criteria (from the VVDS parent survey/catalogues and the addi- 
tional specific criteria considered in choosing our targets), the 
MASSIV sample provides a good representation of "normal" 
star-forming galaxies, with a median SFR ~ 30 Mq yr"', in the 
stellar mass regime 10'' - 10" Mq. 



5. Composite integrated spectra 

5.1. VIMOS spectra 

Composite VIMOS spectra have been built from integrated ID 
spectra of MASSIV galaxies extracted from the VVDS database. 
Each spectrum has been first brought to rest frame (using the 
VVDS redshift) and then normalized to unity over the interval 
2975 - 3025A. The result of the averaging procedure is shown 
in Figure ]7j whichpresents the composite spectrum for the 83 
MASSIV galaxie^(the galaxy VVDS220148046 has been ex- 
cluded due to its redshift z - 2.24) in the rest-frame wavelength 
range ~ 1500 - 4500A. The high S/N of this composite spec- 
trum (~ 30 or typically 5-10 times that of individual spectra 
in the wavelength range 2400 - 3800A) is apparent. Compared 
to individual spectra, we see significantly more absorption lines 
(indicated with blue labels) in the composite as well as faint 
emission lines ([Ne iii]/l3869 and Hy) barely detected in individ- 
ual frames. Composite spectra have also been produced in two 
different stellar mass bins, below (44 galaxies) and above (39 
galaxies) a stellar mass of 1.6 x 10'" Mg, corresponding to the 
median value of the MASSIV sample. These composite spectra 
are shown in Figure ]8] and are used in the next sections to in- 
vestigate any mass dependence on the AGN contamination and 
integrated metallicity. The emission line equivalent widths mea- 
sured in the three composite spectra are listed in Table ]3] 



^ A similar composite spectrum has been built using the 44 MASSIV 
galaxies of the first epoch sample only. No significant difference has 
been found from that for the full sample. 
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Fig. 7. Top: Rest-frame stacked VIMOS spectrum for 
the full MASSIV sample (83 galaxies after exclusion of 
VVDS220148046 at z = 2.24) normalized at the continuum 
level at A - 3000A. The position of main absorption (blue) and 
emission (red) lines is indicated with dotted lines. The middle 
and bottom panels are for the same spectrum but zoomed-in on 
specific wavelength ranges. 



5.2. SINFONI spectra 

Composite SINFONI spectra have been built from integrated ID 
spectra of the "first epoch" sample of MASSIV galaxies. Each 
spectrum has been obtained following the procedure detailed in 
Queyrel et al. (2011), corrected for the "artifical" line broaden- 
ing due to large-scale motions, and brought to rest frame us- 
ing the SINFONI redshift. Spectra are then normalized in order 
to set the Ha peak intensity to unity. The result of the averag- 
ing procedure is shown in Figure [8] which presents the compos- 
ite spectrum for 44 MASSIV galaxie^in the rest-frame wave- 
length range ~ 6450 - 6750A. The high S/N of this compos- 
ite spectrum allows one to see clearly the [Nii]/l/l6548,6584 an 
[S ii]/l/l6717,6731 emission lines which are barely detected in 
a high fraction of individual frames (see Queyrel et al. 2011). 
Composite spectra have been also produced in two different stel- 
lar mass bins, below (23 galaxies) and above (21 galaxies) the 
median value (1.6 x 10'" Mq) of stellar mass. These compos- 
ite spectra are shown in Figure [8] The emission line flux ratios 
measured in the three composite spectra are hsted in Table |3] 

5.3. Physical properties derived from composite spectra 

Thanks to their high S/N, the VIMOS and SINFONI compos- 
ite spectra can be used to measure faint emission lines, such 
as [Nein]^3869, Hy, [Nii]^^6548,6584, or [S ii]/l^6717,6731, 
which are barely detected in a few individual galaxies only. 
These emission lines can then be used to probe some average 



44 galaxies instead of 50 because 4 objects are targeted in 
[Oin]/15007, including the one at ~ 2.24, and 2 galaxies were not de- 
tected 
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Fig. 8. Rest-frame composite VIMOS (left) and SINFONI 
(right) ID spectra, normalized at the continuum level at /I = 
3000A (VIMOS) or at the Ho- intensity level (SINFONI). Top 
panels: spectra are shown for the full (left) and "first epoch" 
(right) MASSIV samples. Spectra corresponding to two stellar 
mass bins are shown in middle (< 1.6 x 10'" M©) and bottom 
(> 1 .6 X 10'" Mq) panels. The position of main absorption (blue) 
and emission (red) lines are indicated with dotted lines. In each 
panel, the number of galaxies used to build the composite spec- 
trum is indicated. 



physical properties of the MASSIV sample, such as the origin of 
the ionizing photons (starburst vs. AGN) and their global metal- 
licity. 



5.3.1. AGN contamination? 

MASSIV targets have been selected to be star-forming galax- 
ies. However, even if type-1 AGNs are clearly identified in the 
VVDS sample thanks to their broad emission lines, the MASSIV 
sample could include a non-negligible fraction of type-2 AGNs 
(typically a few percent). We can use the composite spectra to 
check if low-activity type-2 AGNs are present in our sample, by 
measuring typical emission-line ratios commonly used to distin- 
guish star-forming galaxies from AGNs. 

Using the emission lines measured in the VIMOS compos- 
ite spectrum, namely [On]/13727, [Neiii]/i3869, and Hy (see 
Table [3|, we can use the line ratios based on the [Ne iii]/l3869 
emission line and introduced by Perez-Montero et al. (2007) 
as a first check for a possible contamination by narrow-line 
AGNs. These ratios involve the Wb emission line (see Figure |9]l 
which can be deduced from the Hy emission line assuming a 
Balmer line ratio Hy/H5 - 1.82 typical of star-forming galax- 
ies (Osterbrock 1989). The location in the diagnostic diagrams 
of emission-line ratios measured in the three VIMOS composite 
spectra (i.e. for the full MASSIV sample and for the two stellar 
mass bins) are shown in Figure |9] together with the empirical 
limits (and associated uncertainty regions; see Perez-Montero et 
al. 2007) between star-forming objects and AGNs. Emission-line 
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Table 3. Emission-line equivalent widths and ratios measured in VIMOS and SINFONI composite spectra 



MASSIV Sample 


EW([Oii]^3727) 


EW([Neiii]^3869) 


EW(H7) 


[Nn]^6584/HQ' 


[Sn]ii6717,6731/HQ' 


[Sn]^6717M6731 




[A] 


[A] 


[A] 








Total 


-53.3 ± 2.6 


-5.4 ± 1.1 


-6.7 ± 1.5 


0.154 + 0.012 


0.137 + 0.010 


1.54 + 0.11 


log(M. [Mo]) < 10.2 


-61.2 + 3.1 


-6.8 ± 1.2 


-6.6+1.3 


0.140 + 0.016 


0.112 + 0.015 


1.44 + 0.16 


log(M. [Mo]) > 10.2 


-45.9 ±2.8 


-3.6+1.2 


-6.4+ 1.2 


0.170 + 0.015 


0.178 + 0.016 


1.64 + 0.15 



ratios for the full MASSIV sample and for the two stellar mass 
bins are both very close to the separation limit between star- 
forming galaxies and AGNs. Note, however, that the emission- 
line ratios which are the closest to this limit are those measured 
in the composite spectrum of the low-mass (< 1.6 x 10'" Mq) 
galaxies contrary to what would be expected if the MASSIV 
sample was contaminated by a significant fraction of narrow-line 
AGNs, as active nuclei are preferentially located in the center of 
massive galaxies. The relatively high intensity of [Nenij/13869 
emission-line in low-mass MASSIV galaxies is certainly due to 
their low metallicity (see section 5.3.2 1. MASSIV galaxies pop- 



ulating the high stellar mass range show emission-line ratio typ- 
ical of star-forming galaxies, even taking into account the un- 
certainty region. Their low-intensity [Ne iii]/l3869 emission-line 
can be explained by their higher metallicity (see section 5.3.2 1. 

The SINFONI composite spectra can also be used to 
check for a possible narrow-line AGN contamination in the 
MASSIV sample. Line ratios involving the [Nii]/16584 and 
[S n]/l/l6717,6731 emission lines are listed in Table [5] for the 
three composite spectra. Unfortunately, we do not have ac- 
cess to the [Oiii]/l5OO7/Hj0 line ratios commonly used in stan- 
dard diagnostic diagrams. However, the [N ii]/l6584/Hci' and 
[S ii]/l/l6717,673 l/Hff Hne ratios measured in the MASSIV com- 
posite spectrum (~ 0.15 and 0.14 respectively) are low enough 
to secure the MASSIV galaxies in the star-forming region of the 
diagnostic diagrams (see e.g. Lamareille et al. 2006). The differ- 
ence in the measured line ratios between low-mass and massive 
galaxies is probably due also to a metallicity effect. 

Based on these arguments, we thus conclude for a low con- 
tamination (< 2-3%) by type-2 AGNs, if any, in the MASSIV 
sample. 




log([OII]3727/H(5) 



l°g(02Ne3) 



Fig. 9. Diagnostic diagrams based on [Neinj/13869 emis- 
sion line. In the left-hand panel, the relation between 
the emission-Hne ratios [Ne iii]/l3869/H5 and [Oii]/l3727/H5 
is shown and in the right-hand panel the relation be- 
tween the emission-line ratio [Neiii]/13869/[On]/l3727 and 
the C>2Ne3=([Oii]/13727-i-15.37[Neiii]^3869)/H5 parameter is 
shown. The solid lines represent the analytical division (see 
Perez -Montero et al. 2007) between star-forming galaxies (lower 
left) and AGNs (upper right). The dashed lines represent the 
limit of the bands of uncertainty of 0.15 dex on both sides of 
each relation. Black square, blue point, and red triangle show the 
location of emission-line ratios measured in the composite spec- 
trum of the full MASSIV sample, of galaxies with stellar mass 
lower than 1 .6x10'" Mq, and of galaxies with stellar mass higher 
than 1.6 X 10'" Mq, respectively. Magenta asterisks show the lo- 
cation of individual VVDS galaxies with measured [Ne iii]/l3869 
emission line (Perez-Montero et al. 2009). 



5.3.2. Global metallicity 

An estimate of the average gas-phase oxygen abundance in the 
MASSIV sample can be obtained from the emission-line ratios 
measured in the composite spectra. We refer to Queyrel et al. 
(2011) for a detailed study of spatially-resolved metallicity in 
the individual MASSIV galaxies. 

Using the VIMOS composite spectra, the oxygen abun- 
dance has been estimated from line ratios involving [O ii]/i3727, 
[Nem]/13869, and H6 emission lines together with the calibra- 
tions proposed by Perez-Montero et al. (2007). This method 
has been already used to derive the mass-metallicity relation 
of VVDS galaxies in the redshift range 0.9 < z < 1.24 
(Perez-Montero et al. 2009). We find an oxygen abundance 
12-Hlog(0/H) = 8.47, 8.32, and 8.62 for the full MASSIV sam- 
ple, for galaxies with stellar mass lower than 1.6 x 10'" Mq, and 
for galaxies with stellar mass higher than 1.6 x 10'" Mq, respec- 
tively. 

Some authors (eg. Maier et al. 2006; Liu et al. 2008) claim 
that the typical physical properties (electron density, hardness of 
the radiation field, etc) of the ionized gas in high-redshift star- 
forming galaxies may be different from the ones observed in the 
local universe, which would affect the validity of the metallic- 



ity calibrations based on nearby galaxy samples. Whereas it is 
not possible to check for the ionization conditions through the 
[Oni]/[Oii] line ratio in our MASSIV sample as we are missing 
the [O in]/15007 emission-line, we can make use of the measured 
[S n]6717/[S ii]6731 line ratios listed in Table[3]to infer the aver- 
age electron density (Osterbrock 1989). Taking into account the 
measurement uncertainties, we find electron densities of the or- 
der of ~ 100 cm"^ (independently of the stellar mass range), 
which are typical of nebular conditions in local star-forming 
galaxies. We are thus quite confident in using metallicity cali- 
brations based on local samples. 

Another independent estimate of the global metallicity of 
MASSIV galaxies is obtained from the [N ii]/l6584/Ha' Hne 
ratio measured in the SINFONI composite spectra. The cali- 
bration proposed by Perez-Montero & Contini (2009) is used 
as in Queyrel et al. (2011). We find an oxygen abundance 
12H-log(0/H) = 8.43, 8.39, and 8.46 for the full MASSIV sam- 
ple, for galaxies with stellar mass lower than 1.6 x 10'" Mq, 
and for galaxies with stellar mass higher than 1.6 x 10'" Mq, 
respectively. The global metallicity of MASSIV galaxies esti- 
mated from two different datasets and calibrations are very con- 
sistent and show, as expected from the mass-metallicity relation, 
a higher value for massive galaxies. 



12 



T. Contini et al.: MASSIV. I. Survey description and global properties of the galaxy sample 



6. Comparison with other IFU samples of 
high-redshift galaxies 

For the sake of consistency and to allow for an homogeneous 
and fair comparison between the various galaxy samples de- 
scribed below, we refer to an unique IMF (Salpeter 1955) for the 
derivation of stellar masses and star formation rates. Following 
Bernardi et al. (2010), we thus applied a constant scaling factor 
of +0.25 or +0.15 to those quantities derived using a Chabrier 
(2003; for SINS, OSIRIS and LSD/AMAZE samples) or a "diet" 
Salpeter (for IMAGES sample) IMF in the SED fitting proce- 
dure. 

6.1. IMAGES 

IMAGES (Intermediate MAss Galaxies Evolution Sequence) 
survey targeted galaxies in the Chandra Deep Field South, with 
redshifts z ~ 0.4 - 0.75, Iab < 23.5 and detected [Oii]^3727 
emission lines with rest-frame EW < -15 A (Ravikumar et 
al. 2007). Galaxies were further selected to be of intermediate 
mass, using /-band absolute magnitudes as a proxy for stellar 
mass, such that Mj < -20.3, which roughly corresponds to 
M. > 1.5 X 10'" Mq. The IMAGES galaxies, which are rep- 
resentative of the 7-band luminosity function at these redshifts 
(Yang et al. 2008), have been observed with the multi-integral- 
field unit spectrograph FLAMES/GIRAFFE (with a sampling of 
0.52"/pixel) at the VLT 

The final IMAGES sample contains a total of 63 galaxies 
with IFU resolved kinematics. The redshift range of these galax- 
ies extends from z ~ 0.4 to z ~ 0.75 with a median value of 



0.61 (see Figs. 11 and 10 1. The stellar masses have been es- 



timated for 40 IMAGES galaxies (Puech et al. 2010; Neichel 
et al. 2008) using a "diet" Salpeter IMF and Bell et al. (2003) 
simplified recipes for deriving stellar mass from 7-band lumi- 
nosity. The distribution of stellar masses for the IMAGES sam- 
ple is shown in Fig 
to 2 X 10" M, 
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Stellar masses range from ~ 6.3 x 10^ 
10 with a median value of 3.5 x 10"^ M©. There 
is no estimate of the SED-based SFR for the IMAGES sample. 
However, the "total" SFR has been computed (see Puech et al. 
2007 for details) as the sum of SFRuv, derived from the 2800A 
luminosity, and SFRir, derived from Spitzer/MIPS photometry 
at 24/im using the Chary & Elbaz (2001) calibration between 
rest-frame 15//m flux and the total IR luminosity. The calibra- 
tions of Kennicutt (1998), which rely on a Salpeter IMF, have 
been used to convert both UV and IR luminosities into SFRs. 
The distribution of total SFR for 36 IMAGES galaxies is shown 



in Fig. 13 Total SFR ranges from ~ 5 to 80 Mg yr , with a 



median value of 10 Mq yr . 
6.2. SINS 

With a total of 62 galaxies at 1.4 < z < 2.6 detected in Ha, the 
SINS (Spectroscopic Imaging survey in the Near-infrared with 
SINFONI; Forster Schreiber et al. 2009) survey, carried out with 
SINFONI at the VLT, is one of the largest surveys of near-IR 
integral field spectroscopy to date. The SINS galaxies were se- 
lected from the spectroscopically-confirmed subsets of various 
imaging surveys in the optical, near-IR, and mid-IR. The pho- 
tometric selection of the parent samples encompassed a range 
of star-forming populations at high redshift, including optically- 
selected "BX/BM" galaxies at z ~ 2 - 3, and infrared-selected 
galaxies at z ~ 1.5 - 2.5, with a majority of "BzK" objects. The 
SINS targets were further selected to have an integrated emission 
line (Ha or [O iii]/l5007 in a few cases) flux of at least 5 x 10"'^ 
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Fig. 10. Redshift distribution of the MASSIV sample (top panel) 
compared with other samples of high-z galaxies observed with 
IFUs: the IMAGES sample (2nd panel from the top), the SINS 
Ha sample (3rd panel from the top), the OSIRIS sample (4th 
panel from the top), and the LSD/AMAZE sample (bottom 
panel). 



ergs s"' cm"^. Although with some bias towards the bluer part of 
the galaxy population compared to purely ^T-band selected sam- 
ples at similar redshifts, the SINS Ha sample provides a reason- 
able representation of massive actively star-forming galaxies at 
z~2. 



The redshift histogram of the SINS galaxies (see Fig. 10 1 
shows a bimodal distribution, peaking at z ~ 1.6 and z ~ 2.3, 
and has a median value of z ~ 2.17. The reason for the bi- 
modal distribution is simply the result of the requirement that 
the emission line of interest (primarily Ha) falls within either 
the H or K band atmospheric windows. The stellar masses have 
been estimated for 59 SINS galaxies using a standard SED fitting 
procedure based on extensive photometry (see details in Forster 
Schreiber et al. 2009), and Bruzual & Chariot (2003) stellar pop- 
ulation synthesis models, with a constant SFR, solar metallicity, 
and a Chabrier (2003) IMF. The distribution of stellar masses 



for the SINS sample is shown in Figure 12 Stellar masses range 
from ~ 3 X 10** to 5 X 10" Mq with a median value of 4.6 x 10'" 
Mg. SED-based star formation rates have also been derived from 
stellar population models for this SINS sub-sample of 59 galax- 



ies. Their distribution is shown in Figure 13 SED-based SFR 
ranges from ~ 1 to 1300 Mq yr"', with a median value of 129 

Moyr-'. 

6.3. OSIRIS sample 

Galaxies in the so-called "OSIRIS sample" (Law et al. 2007, 
2009; Wright et al. 2007, 2009) have been selected from the rest- 
frame UV color-selected catalogue of Steidel et al. (2004), and 
are all confirmed with optical spectroscopy using LRIS on the 
Keck telescope. The Steidel et al. survey primarily focused on 
z ~ 2 - 3 galaxies, but has also used the two-color technique 
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Sample 




Redshift 




log(Stellar Mass) [Mq] 


log(SFRsED)' 


[Mo yr-'] 




Range 


Median ± a 


N 


Median ± cr 


A' 


Median ± cr 


N 


MASSIV 


0.9 - 1.8 


1.33 ± 0.13 


83 


10.15 ±0.30" 


83 


1.50 ±0.31 


83 










10.21 ±0.27^ 


83 






IMAGES 


0.4-0.75 


0.62 ± 0.07 


63 


10.54 ± 0.24 


40 


0.98 ± 0.29 


36 


SINS 


1.4-2.6 


2.17 ±0.36 


62 


10.66 ±0.37 


59 


2.11 ±0.31 


59 




1.4-2.0 






10.76 ±0.38-* 


18 


2.43 ± 0.48"* 


18 


OSIRIS 


1.5-3.3 


2.17 + 0.32 


20 


10.29 ± 0.43 


20 


1.92 ±0.35 


14 




1.5-2.0 






10.20 ±0.41' 


7 






LSD/AMAZE 


2.6-3.8 


3.29 ±0.18 


18 


10.31 ±0.30 


18 


2.25 ±0.16 


9 



' Except for the IMAGES sample, see text for details 

^ Stellar mass derived from SED fitting using the PDF method 

^ Stellar mass derived from SED fitting using the "best fit" method 

Stellar mass and SFR derived from SED fitting for the SINS z < 2 subsample only 
' Stellar mass derived from SED fitting for the OSIRIS z <2 subsample only 




6 2 
Look-back time [Gyr] 

Fig. 11. Evolution of the cosmic star formation rate density as a 
function of look-back time and redshift (adapted from Hopkins 
2006). The redshift range of MASSIV (0.9 <z < 1.8, magenta 
box) is compared with other major IFU surveys of distant galax- 
ies: IMAGES (z ~ 0.4-0.75), SINS/OSIRIS (z ~ 1.4-2.6), and 
LSD/ AMAZE (z ~ 2.6 - 3.8) . The relative height of each boxes 
is proportional to the samples size. 



introduced in Adelberger et al. (2004) to detect lower redshift 
z ~ 1.6 star-forming counterparts within the traditional redshift 
desert (Steidel et al. 2004). The heterogeneous OSIRIS sample 
has been built following a variety of criteria. Some targets were 
deliberately selected for their young stellar population and cor- 
respondingly small stellar masses, some for their old ages and 
large stellar masses, and some for other reasons including com- 
plex or multi-component rest-frame UV morphologies (Law et 
al. 2007), strong detections in Ha narrowband surveys, unusual 
spectral features, or previous acquisition of long-slit kinematic 
data. However, given the relatively shallow OSIRIS /T-band sen- 
sitivity, the most common criterion applied was based on the 
emission line flux (> 5 x 10"'^ ergs s"' cm"^) previously mea- 
sured with long-slit spectroscopy (Erb et al. 2006). IFU obser- 
vations have been obtained at the Keck telescope using NIR 
integral-field spectrograph OSIRIS with the AO/LGS system. 




10 11 12 

log(Stellar Mass) [Msun] 

Fig. 12. Stellar mass distribution of the MASSIV sample (top 
panel). The filled (empty) histogram is for PDF ("best fit") 
masses. This distribution is compared with other samples of 
high-z galaxies observed with IFUs: the IMAGES sample (2nd 
panel from the top), the SINS Ha sample (3rd panel from the 
top), the OSIRIS sample (4th panel from the top), and the 
LSD/AMAZE sample (bottom panel). The filled histograms in 
the 3rd and 4th panel are restricted to z < 2 galaxies. For 
comparison between the different samples, a scaling factor of 
-hO.25 dex (SINS, OSIRIS, and LSD/AMAZE samples) or -h0.15 
dex (IMAGES sample) has been applied to the published stellar 
masses. 



The OSIRIS sample contains a total of 20 galaxies with Ha 
measurements. The redshift range of these galaxies extends from 
1.5 to z ~ 3.3 (see Figs.fTT|and 



z ~ 1.3 to z ~ 3.3 (see Figs 
bution peaking at z ~ 1.6 and z 



10 1, with a bimodal distri- 
2.2, and a median value of 



z ~ 2.17. The reason for the bimodal distribution is the same as 
for the SINS sample: the emission line of interest (primarily Ha) 
has to fall within either the H or K band atmospheric windows. 
The stellar masses have been estimated using a standard SED 



14 



T. Contini et al.: MASSIV. I. Survey description and global properties of the galaxy sample 



T — 1 — I — 1 — 1 — 
MASSIV sample 



H h 



H 1 h 



H 1- 



H 1 1 — ^ 



-1 1 1 1 

IMAGES sample 



H 1- 




V—^ 1 H 



H 1 1 1 1 1 h 



SINS Ha sample 



H 1 1 1 1 1 1- 

: OSIRIS sample 




H 1 1 1 r-p l i l—T^ 



H 1- 



: LSD/AMAZE sample 



R=P — I l-C^p — I H 



12 3 

log(SFRgE„) [Msun/yr] 

Fig. 13. SED-derived SFR distribution of the MASSIV sample 
(top panel) compared with other samples of high-z galaxies ob- 
served with IFUs: the IMAGES sample (2nd panel from the top), 
the SINS Ha sample (3rd panel from the top), the OSIRIS sam- 
ple (4th panel from the top), and the LSD/ AMAZE sample (bot- 
tom panel). The filled histograms are restricted to z < 2 galax- 
ies. For comparison between the different samples, a scaling fac- 
tor of +0.25 dex (SINS, OSIRIS, and LSD/AMAZE samples) or 
+0.15 dex (IMAGES sample) has been applied to the published 
star formation rates. 



fitting procedure based on optical/NIR (UGRJK) and MIR/FIR 
(Spitzer IRAC and MIPS) photometry, and Bruzual & Chariot 
(2003) stellar population synthesis models, with a constant SFR, 
solar metallicity, and a Chabrier (2003) IMF. The distribution 



of stellar masses for the OSIRIS sample is shown in Figure 12 
Stellar masses range from ~ 2 x 10^ to 10" M© with a median 
value of 2.0 X 10'" Mq. SED-based star formation rates have 
also been derived from stellar population models for a fraction 
(14 out of 20 galaxies) of the OSIRIS sample . Their distribution 
is shown in Figure 13 SED-based SFR ranges from ~ 8 to 700 



Mq yr with a median value of 83 Mq yr 



6.4. AMAZE & LSD 

AMAZE & LSD are two surveys conducted with SINFONI, 
with the main aim of studying the metallicity and dynamics 
of high-redshift (z > 3) galaxies. The AMAZE (Assessing the 
Mass-Abundance redshift Evolution; Maiolino et al. 2008) sam- 
ple consists of 30 Lyman Break Galaxies in the redshift range 
3 < z < 4.8 (most of which at z ~ 3), with deep Spitzer/IRAC 
photometry (3.6 - 8 fim). These galaxies were observed with 
SINFONI in seeing-limited mode. LSD (Lyman-break galaxies 
Stellar populations and Dynamics; Mannucci et al. 2009) is a 
smaller sample of 10 LBGs at z ~ 3 with available Spitzer 
and HST imaging. For LSD, SINFONI observations were per- 
formed with the aid of adaptive optics in order to improve spa- 
tial resolution since this project was aimed to obtain spatially- 
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Fig. 14. SED-derived star formation rate as a function of stel- 
lar mass. The MASSIV sample (black filled squares) is com- 
pared with other major IFU surveys of distant galaxies: IMAGES 
(z ~ 0.4-0.75, blue circles), SINS (z ~ 1.4-2.6, magenta trian- 
gles), OSIRIS (z~ 1.5 - 3.3, green crosses), and LSD/ AMAZE 
(z ~ 2.6 - 3.8, red asterisks). The lines represent the empiri- 
cal relations between SFR and stellar mass for different redshifts 
between z = and z = 3 following the analytical expression 
given in Bouche et al. (2010). For comparison between the dif- 
ferent samples, a scaling factor of +0.25 dex (SINS, OSIRIS, 
and LSD/AMAZE samples) or +0.15 dex (IMAGES sample) has 
been applied to the published stellar masses and star formation 
rates. 



resolved spectra for measuring kinematics and gradients in emis- 
sion lines. 

The combined LSD/AMAZE sample (Gnerucci et al. 2011) 
contains a total of 35 galaxies with [Oiii]/I5007 emission line 
measurements, except for the nearest z ~ 2.6 galaxy for which 
Ha has been targeted. The redshift range of these galaxies ex- 
tend from z ~ 2.6 to z ~ 3.8 with two galaxies only at z < 3 (see 
Fig. 
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and a median value of ~ 3.3. The stellar masses have 
been estimated using a standard SED fitting procedure based 
on optical/NIR (UGRJK) and Spitzer IRAC photometry, and 
Bruzual & Chariot (2003) stellar population synthesis models, 
with a smooth exponentially decreasing SFR, a Z = 0.2 Zq 
metallicity, and a Chabrier (2003) IMF. The distribution of stel- 
lar masses for the LSD/AMAZE sample is shown in Figure 12 
Stellar masses range from ~ 2 x 10^ to 10'^ Mq with a me- 
dian value of 2.0 x 10'" Mq. SED-based star formation rates 
have also been derived from stellar population models for the 
AMAZE sample only (9 galaxies) . Their distribution is shown 
in Figure 13 SED-based SFR ranges from ~ 60 to 750 Mq yr"'. 



with a median value of 178 Mq yr 



6.5. Comparison 



Table |4] summarizes the redshift range probed by the different 
IFU surveys presented in the previous sections, together with 
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the median values of redshift, stellar mass and SED-based SFR, 
allowing a direct comparison with MASSIV. 

The redshift range probed with the MASSIV sample is com- 



pared with the other IFU samples in Figures 10 and 1 1 MASSIV 
covers a redshift range (z ~ 0.9 - 1.8) which is intermedi- 
ate between the high-z surveys (z ~ 2-3, SINS, OSIRIS 
and LSD/AMAZE) and the IMAGES sample (z ~ 0.4 - 0.8). 
However, the MASSIV, OSIRIS and SINS surveys are probing a 
common redshift range at z ~ 1.3-1.8, whereas LSD/AMAZE 
targets only galaxies at the highest redshifts (z ~ 2.6 - 3.8). 
MASSIV could in principle follow the evolution of galaxy kine- 
matics and mass assembly at a crucial epoch corresponding to a 
transition phase between non-relaxed systems observed at high 
redshifts and more evolved and stable galaxies which populate 
the Hubble sequence at low redshifts. 

The stellar mass distribution of the MASSIV galaxies is 
compared with the other IFU samples in Figure 12 The adopted 
stellar masses for the MASSIV sample are those derived with the 
PDF method (see Section [43]). However, in order to be able to 
compare with other surveys, we also derived the stellar masses 
with the so-called "best fit" (BE hereafter) method which is the 
one adopted for the other IFU samples. Both stellar mass dis- 
tributions (PDF and BE) are thus shown in Figure 12 The stel- 
lar mass range probed by MASSIV is rather similar to the one 
covered by the IMAGES, SINS, OSIRIS, and LSD/AMAZE sur- 
veys, extending from a ~ 10^ to 10" Mq. The median value of 
the MASSIV sample (1.6 x 10'° Mq with the BE method) is sim- 
ilar to the median values of the IMAGES/OSIRIS/LSD/AMAZE 
samples (~ 2 x 10'" Mq) but lower by a factor of ~ 3 (at a con- 
fidence level > 3cr) than the median value of the SINS sample 
(4.6 X 10'° Mo). 

The SED-based SFR distribution of the MASSIV sample is 
also compared with other IFU samples in Figure 13 The lowest 



values of SFR probed by MASSIV (a few Mq yr^^Ts rather sim- 
ilar to the ones targeted by the other surveys, except for the high- 
est redshift surveys LSD/AMAZE starting at ~ 60 Mq yr"' . The 
median value of the MASSIV sample (32 Mq yr"') is, however, 
significantly lower (at a confidence level > 3cr) than the median 
values of the OSIRIS (83 Mq yr"'), SINS (129 Mq yr"'), and 
LSD/AMAZE (178 Mq yr"') samples. 

In Figure[T4]we compare the location of the MASSIV galax- 
ies to the other IFU samples in the SFR vs. stellar mass rela- 
tionship. As already stated in sect. 4.2 the MASSIV sample is 
representative of the overall population of star-forming galaxies 
at z ~ 1-2. For a given stellar mass, the level of star forma- 
tion probed by MASSIV is globaly lower (at a confidence level 
> 3cr) that the one of SINS and LSD/AMAZE targets which 
reach easily SFRs higher than ~ 500 Mq yr '. On the opposite, 
the IMAGES sample shows much lower SFRs (< 100 Mq yr ') 
for a given stellar mass as expected by the star formation "main 
sequence" defined at z ~ 0.4 - 0.8 (see Bouche et al. 2010). 

Compared to other high-z IFU surveys, the main advantage 
of MASSIV is its representativeness for star-forming galaxies 
with SFR > 5 Mq yr '. Indeed, MASSIV has been drawn from 
VVDS, the most unbiased and representative sample of high- 
z star-forming galaxies available so far. On the contrary, sur- 
veys like SINS, OSIRIS and LSD/AMAZE targeted galaxies se- 
lected in various (and hence heteregeneous) color-selected sam- 
ples (LBG, BzK, sub-mm galaxies, etc), thus probing mainly 
galaxies at a relatively high star formation level (see Table|4|. 

One could argue that this difference in SFR is largely at- 
tributed to a combination of higher redshift ranges covered by 
other IFU surveys (except IMAGES) and the limits imposed on 
the minimum (expected) line flux or SFR to ensure detection. 



But if we restrict the SINS sample to the 1.4 - 2.0 redshift range 
in common with MASSIV, the SINS galaxies have still a median 
SFR (~ 270 Mq yr"'; 18 objects) significantly higher than the 
MASSIV one (~ 60 Mq yr"'; 27 objects). The same is true for 
the expected lower limit on SFR imposed by a selection based 
on Hne flux ([Oii]^3727 for MASSIV and Ho- for a fraction of 
SINS galaxies). Using the quoted lower limit on Ho- flux used 
for the selection of SINS galaxies (= 5 x 10"'^ ergs s"' cm"^; 
Forster Schreiber et al. 2009), we end up with a minimum SFR 
of 12 Mq yr"', similar to the value derived for MASSIV (8 Mq 
yr"' for [Oii]/l3727 line flux = 1.5 x 10"'^ ergs s"' cm"^), us- 
ing the Argence & Lamareille (2009) consistent calibrations to 
derive the SFR from line luminosities. 

The total number of 84 galaxies included in MASSIV is 
also a clear advantage compared to surveys like OSIRIS and 
LSD/AMAZE which are limited to ~ 20 galaxies each. The 
MASSIV sample indeed allows us to probe different mass and 
SFR ranges while keeping enough statistics in each category 
in order to derive robust conclusions on galaxy kinematics and 
mass assembly. 

7. Conclusions 

This paper presents the basic physical properties of the full 
MASSIV sample containing 84 star-forming galaxies drawn 
from the VVDS in the redshift range 0.9 < z < 2.2. SINFONI 
observations are now completed and the analysis of a first sam- 
ple of 50 MASSIV galaxies is presented in related papers fo- 
cusing on the kinematical classification (Epinat et al. 2011), 
the evolution of fundamental scaling relations such as the bary- 
onic Tully-Fischer one (Vergani et al. 2011), and the spatially- 
resolved metallicity (Queyrel et al. 2011). The selection criteria 
are mainly based on [Oii]/l3727 emission line strength, mak- 
ing the MASSIV sample a fair representation of "normal" star- 
forming (with a median SFR ~ 30 Mq yr"') galaxies in the 
stellar mass regime 10'' - 10" Mq. Analysis based on com- 
posite VIMOS and SINFONI spectra reveals that the contami- 
nation by type-2 AGNs, if any, is very low (< 2-3%) and that 
the integrated metallicity of the galaxies follow the well-known 
mass-metallicity relation. Compared to other existing high-z IFU 
surveys, we conclude that the main advantages of MASSIV are 
its representativeness (flux-selected from the magnitude-selected 
VVDS sample), and sample size allowing us to probe different 
mass and SFR ranges while keeping enough statistics in each 
category in order to derive robust conclusions on galaxy kine- 
matics and mass assembly. 
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Table 5. Global properties of the MASSIV sample. Column 1 lists the VVDS identification number. Columns 2 and 3 give the RA and DEC 
coordinates. Columns 4 and 5 list the galaxy redshift and associated flag (see Le Fevre et al. 2005) as measured in VIMOS spectra. Column 6 
indicates the depth (D=Deep, W=Wide, UD=Ultra-Deep) of the VVDS parent sample. Column 7 gives the galaxy stellar mass and associated 
uncertainty (in log of solar masses). Colurrm 8 lists the galaxy SED-based star formation rate and associated error (in log of Mq yr '). Columns 
9 and 10 indicate the [On]/l3727 equivalent width (in A) and flux (in 10"'^ ergs s"' cm"^) and associated uncertainties as measured in VIMOS 
spectra 
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+ 


0, 
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.80 ±0.41 
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Table 5. continued. 



VVDS ID 


RA(2000) 


Dec(2000) 


Redshift 


Zflag 


Sample 




SFRsED 


EW[oii] 


F[Oii] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


910184233 


02:27:18.779 


-4:33:49.288 


1.4614 


2 


UD 


9.55 ± 0.35 


QO 


-54.0 ± 4.2 


2.1 ± 0.4 


910186191 


02:26:17.891 


-4:33:34.027 


1.5399 


23 


UD 


11.27 ±0.45 


9 99+u.ou 

^•^^-0.32 






910187744 


02:27:17.235 


-4:33:20.052 


1.3191 


3 


UD 


10.07 ± 0.20 


1.44 ±0.18 


-79.7 ± 7.2 


5.6 ± 0.5 


910191357 


02:25:26.617 


-4:32:51.025 


1.2863 


3 


UD 


9.99 ± 0.22 


1.13 ±0.39 


-48.0 ± 4.5 


1.3 ±0.2 


910193711 


02:25:46.278 


-4:32:34.066 


1.5523 


4 


UD 


9.99 ± 0.18 


1.95 ±0.31 






910195040 


02:27:32.118 


-4:32:24.490 


1.7712 


4 


UD 


10.58 ± 0.20 








910207502 


02:26:17.248 


-4:30:47.592 


1.4656 


4 


UD 


10.31 ±0.33 


1.97!«S 


-55.0 ±5.1 


4.0 ± 0.4 


910224801 


02:25:37.690 


-4:28:32.606 


1.3201 


3 


UD 


9.96 ± 0.20 


1.34 ±0.50 


-56.3 ± 4.8 


3.3 ± 0.3 


910232719 


02:26:42.439 


-4:27:26.554 


1.7592 


3 


UD 


10.37 ± 0.22 


1.96 ±0.32 






910238285 


02:27:16.371 


-4:26:44.682 


1.5405 


3 


UD 


10.92 ±0.13 


1.89 ±0.18 






910247797 


02:27:11.027 


-4:25:31.559 


1.5267 


3 


UD 


9.77 ± 0.35 


1 30+0.53 

1.55 ± 0.40 






910250031 


02:26:33.021 


-4:25:08.908 


1.6609 


2 


UD 


10.15 ± 0.20 






910254325 


02:27:29.241 


-4:24:40.136 


1.2818 


2 


UD 


10.22 ± 0.27 


1.36 ± 0.41 


-38.8 ± 3.9 


1.9 ±0.2 


910259245 


02:26:44.889 


-4:23:56.486 


1.5174 


3 


UD 


10.79 ±0.12 


1 Q'7+U.Uu 






910261247 


02:26:52.706 


-4:23:45.074 


1.4262 


4 


UD 


10.57 ± 0.14 


1.70 ±0.39 


-70.0 ± 6.5 


4.0 ± 0.3 


910262816 


02:27:13.437 


-4:23:29.756 


1.4807 


3 


UD 


11.26 ± 0.15 


1 50 
1.50 ± 0.36 


-65.4 ±5.8 


1.6 ± 0.2 


910266034 


02:27:09.381 


-4:23:04.596 


1.5642 


4 


UD 


10.07 ± 0.20 






910274060 


02:25:48.031 


-4:21:58.842 


1.5680 


3 


UD 


10.09 ± 0.21 


1 09+0.55 

^■■'^-0.39 






910276733 


02:26:24.178 


-4:21:42.091 


1.3390 


4 


UD 


10.14 ±0.19 


1.28 ±0.55 


-66.6 ± 4.8 


5.1 ±0.5 


910279515 


02:25:36.246 


-4:21:15.772 


1.3983 


3 


UD 


10.79 ±0.14 


1.90 ±0.07 






910279755 


02:25:48.722 


-4:21:15.610 


1.3127 


4 


UD 


9.96 ± 0.24 


1 34+0.45 


-56.9 ± 5.7 


5.8 ±0.5 


910286831 


02:26:17.403 


-4:20:15.061 


1.4696 


3 


UD 


10.28 ± 0.27 


^ ~ 1 A AC 


-70.0 ± 6.0 


3.8 ±0.3 


910296626 


02:25:50.696 


-4:19:00.556 


1.3570 


2 


UD 


9.87 ± 0.28 


1.24 ±0.53 








AT. 97. 97 /too 


A. 1 C-'J'J 1 9/1 


i .0 / Vi 


J 


T Tr> 




1 99 -1- n ^9 






910337228 


02:26:50.445 


-4:13:19.920 


1.3972 


3 


UD 


9.96 ± 0.20 


1.34 ±0.40 


-107.4 ±8.2 


6.2 ± 0.5 


910340496 


02:27:14.061 


-4:12:54.457 


1.3981 


3 


UD 


10.03 ± 0.25 


1.20 ±0.30 


-67.0 ± 5.6 


3.2 ± 0.3 


910360676 


02:27:10.203 


-4:10:25.082 


1.7230 


4 


UD 


10.56 ±0.14 


2 37+0.18 






910370574 


02:26:55.915 


-4:09:02.995 


1.6632 


2 


UD 


10.09 ± 0.24 






910371309 


02:27:12.473 


-4:08:59.158 


1.7926 


3 


UD 


10.44 ± 0.30 


2 J 3 +8:^? 






910377628 


02:27:18.456 


-4:08:06.209 


1.4833 


2 


UD 


10.28 ± 0.30 


1 co+8:A§ 

^••"*-0.15 


-110.7 ±8.7 


2.1 ±0.2 



For this galaxy, we indicate the redshift derived from Ha in the SINFONI data (see text for details) 



